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Abstract 

We present a derivation of the low-energy effective meson La- 
grangian of the Nambu - Jona-Lasinio (NJL) model on the basis of 
Schwinger's proper-time regularization of the one-loop fermion deter- 
minant. We consider the case in which the SU{2) x SU (2) chiral 
symmetry of the NJL Lagrangian is broken by the current quark mass 
matrix with rh u ^ rhd- The non-degeneracy of d and u masses de- 
stroys one of the most crucial features of the proper-time expansion - 
the chiral-invariant structure of Seeley - DeWitt coefficients. We show 
however that systematic resummations inside the proper-time expan- 
sion are still possible and derive a result which is in full agreement 
with the chiral Ward - Takahashi identities. 

PACS number(s): 12.39.Fe, 11.30.Rd. 
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1 Introduction 



In recent papers |T], we have shown how to determine systematically the 
low-energy structure of the Nambu - Jona-Lasinio (NJL) model |§ on the 
basis of Schwinger's proper time regularization of the one-loop fermion de- 
terminant 0, [5]. We have considered the cases with linear and non-linear 
realizations of explicitly broken SU{2) x SU{2) chiral symmetry. We have 
shown that in the presence of the explicit chiral symmetry breaking term in 
the Lagrangian, the standard definition of In | det D\ in terms of a proper-time 
integral 

In | det D\ = -\ jT ^p(T, A 2 )Tr (e~ TDtD ) (1) 

modifies the explicit chiral symmetry breaking pattern of the original quark 
Lagrangian and needs to be corrected in order to lead to the fermion deter- 
minant whose transformation properties exactly comply with the symmetry 
content of the basic Lagrangian. The reason for the necessary modification 
of the result obtained on the basis of formula ([I]) is directly connected with 
the restrictions imposed by the chiral Ward - Takahashi (WT) identities. 

In the present paper we extend this framework to the case in which the 
SU(2) x SU(2) chiral symmetry of the NJL Lagrangian is broken by the 
current quark mass matrix with m u ^ m^. The non-degeneracy of d and 
u masses introduces a rather large violation of isotopic spin conservation 
and destroys one of the most crucial features of the proper time expansion 
- the chiral-invariant structure of Seeley - DeWitt coefficients. It means 
that additionally to the already known problem of correcting the standard 
definition ([]]) by the functional P being proportional to the current quark 
masses |l| one has to find a method which replaces the standard asymptotic 
expansion of the heat kernel in terms of Seeley - DeWitt coefficients by a 
new expansion in terms of chiral-invariant structures. As we shall show, 
systematic resummations inside the proper-time expansion are possible. We 
describe here how calculations can be organized in this, at the first sight 
quite fuzzy, situation. The resulting series for the heat kernel is not anymore 
a proper-time expansion. We calculate its first terms and present the method 
to determine higher orders. 

The spontaneous breakdown of the global SU(2) x SU(2) chiral symmetry 
in the NJL model is a consequence of the fact that the Schwinger - Dyson 
equations for the fermion propagators ("gap" -equations) have nontrivial so- 
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lutions with m u ^ rrid ^ 0. After bosonization the same equations occur as 
conditions for the minimum of the effective potential. This generates a corre- 
sponding redefinition of the scalar fields, from which it follows that the part 
of the effective Lagrangian which is linear in the scalar fields vanishes (due to 
the gap-equations), for excitations about the real vacuum state. As we shall 
show, this immediately affects the isospin symmetry breaking part of the La- 
grangian and breaks the chiral WT identities at this level. As a consequence 
we need to invoke new counterterms which would have not been required if 
there had not been spontaneous symmetry breakdown. These counterterms 
have a simple structure and are completely fixed by chiral symmetry. 

This article presents a new systematic expansion procedure for the heat 
kernel of the one loop fermion determinant for the case m u ^ rhd. It differs 
from the methods which have already been used in the literature in the 
same context, for example in papers |7|. The result is also different. In 
Sec. II we discuss the Lagrangian of the NJL model and show that chiral 
SU(2) x SU(2) transformations of quark fields dictate the transformation 
laws of the auxiliary bosonic fields. These collective variables are necessary 
to rearrange the four-quark Lagrangian of the NJL model in an equivalent 
Lagrangian which is only quadratic in the quark fields. The chiral symmetry 
is broken explicitly by the current quark mass matrix with fh u ^ rhd- We 
use the classical equations of motion for the pseudoscalar and scalar fields 
to rewrite the variation of the NJL Lagrangian with respect to the action of 
the chiral group in terms of meson fields. In Sec. Ill we discuss in detail 
the chiral WT identities for the case under consideration. We show how 
one can integrate these identities to get the symmetry breaking part of the 
effective Lagrangian without explicit calculations of the quark determinant. 
In Sec. IV we calculate the fermion determinant. We show how to define 
it for the case in which explicit symmetry breaking takes place and isospin 
symmetry is broken. We calculate the first terms of the new expansion of the 
heat kernel in full detail. We derive the corresponding correcting polynomial 
from the functional P and show that it is completely fixed by the symmetry 
requirements. The effective meson Lagrangian is obtained at the end of this 
section. In Sec. V we introduce new variables for pseudoscalar and scalar 
fields in order to describe the physical meson states and to get the meson 
mass spectrum. The concluding remarks are given in Sec. VI. Finally we 
show in the Appendix some technical details in our treatment of the heat 
kernel exponent with different masses for u and d quarks. 
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2 From quarks to mesons: the infinitesimal 
chiral transformations 



Consider the effective quark Lagrangian of strong interactions with only light 
u and d quark fields which is invariant under a global colour SU (N c ) sym- 
metry 

G 

C = g^Vfy - rh)q + -\{qT a qf + {qi^T a qf). (2) 

Here q is a flavor doublet of Dirac spinors for quark fields q = (u, d). The 
fermion field carries both flavor and color indices. Summation over the color 
indices is implicit. For the tau 2x2 matrices, r a , where a = 0, 1, 2, 3, and 
tr(T a r&) = 25 a b, we shall use the following notation: To = 1 and 7*, (i = 1, 2, 3) 
are the standard Pauli matrices. The manifest chiral symmetry breaking oc- 
curs via the current quark mass matrix: rh = diag(m u , m^) where m u ^ m<j- 
Without this term the Lagrangian (|2|) would be invariant under global chiral 
SU (2) x SU (2) symmetry. We restrict our consideration to the case of only 
scalar and pseudoscalar four-quark interactions with the coupling constant 
G being strong enough to create a stable vacuum state with a nontrivial so- 
lution of the gap-equations corresponding to spontaneous breakdown of the 
chiral 577(2) x SU(2) symmetry. 

The transformation law for the quark fields is the following 

Sq = i(a + 'j 5 P)q, 5q = -iqia-^p), (3) 

where parameters of global infinitesimal SU(2) x SU(2) chiral transforma- 
tions are chosen as a = a^, (3 = (3^. Therefore the Lagrangian C trans- 
forms according to the law 

5C = iq([a,m}-'j 5 {P,m})q. (4) 

It is clear that nothing must destroy this property of the model (we are not 
considering anomalies here), otherwise the content of the theory would be 
changed. In particular it means that the symmetry breaking pattern of the 
bosonized NJL Lagrangian should be the very same, being written only in 
terms of collective fields. 

To make the last statement more rigorous let us consider the path in- 
tegral representation for the generating functional Z in a theory with the 
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Lagrangian density (0). To integrate over anticommuting c- number quark 
fields in Z one has to introduce color singlet collective bosonic variables in 
such a way that the action becomes bilinear in the quark fields and the 
integration becomes trivial 



J VqVqVs a Vp a exp ji J d 4 x 



(5) 



We suppress external sources in the generating functional Z and assume 
summation over repeated flavor indices. In order not to destroy the symmetry 
of the basic quark Lagrangian £ one has to require from the new collective 
variables that 

5(sl+pl) = (6) 

Since there are no kinetic terms for s a and p a , these fields are auxiliary. Let 
us replace the variables in Z: s a ,p a — > cr a , 7r , 

s a = a a -m a + G(qT a q), (7) 

p a = n a - G(qij 5 T a q). (8) 

Here we put 

m u + rh d . rh u -rh d 
m = , rrii = b i3 . (9) 

Requirement (|B|) together with the chiral transformation laws for the quark 
fields (^0 lead to the transformation laws for the new collective fields: 

5a = i[a, a — rh] — {/3, 7r}, 5tt = i[a, ir] + {(3, a — rh}, (10) 

where 7r = n a T a , a = a a r a , and rh = rh a r a . In this way the transformation 
laws of the quark fields with respect to the action of chiral 577(2) x SU(2) 
group define the transformation laws of the collective bosonic fields. However 
Eq. (|T0"|) is not yet the final form which we need. In the N JL model the vacuum 
is not invariant under chiral SU(2) x SU(2) transformations. To explore the 
properties of the spontaneously broken theory one has to define new scalar 
fields with vanishing vacuum expectation values, i.e., one has to rewrite the 
Lagrangian of the theory in terms of shifted fields: 

O — ► (T -f- 771, (11) 
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where 

m u + m d m u -m d 
m = m a T a , m = , rrn = b i3 , (12) 

are the masses of constituent quarks. The chiral transformation properties 
of scalar and pseudoscalar fields change correspondingly to the new ones: 

5a = i[a,a + A] - {/?, tt}, 5tt = i[a, vr] + {(3, a + A}, (13) 

where A = m — rh. This is the final form of infinitesimal chiral transforma- 
tions for collective mesonic excitations around the nontrivial vacuum state 
in the NJL model. 

In accordance with replacements (|7|) and (||) we obtain the mixed form 
for the NJL Lagrangian which includes both quarks and collective degrees of 
freedom 

C(q, q, a a , vr a ) = qDq - —[{a a + A a ) 2 + 7c 2 a ], (14) 

where 

D = ij^d^ — m — o + 275^. (15) 

The Euler - Lagrange equations for mesonic fields take the form of constraints 

7r a = Gqi-f 5 r a q, a a = -A a - Gqr a q. (16) 

By using them one can rewrite Eq.@ in terms of collective degrees of free- 
dom: 

5C = ^Sa a . (17) 

One concludes that the symmetry breaking pattern of the NJL model is not 
affected by the occurence of spontaneous symmetry breakdown. 

The derivation of the action for collective fields 7r a and o a in the NJL 
model is reduced (after integrating out the quark fields) to the calculation of 
the functional determinant of the operator D 

5 coll = / dixC C o\\ = ~ nndetD -^GJ d * x [ (cTa + Aa)2 + A ■ (18) 

This is a great advantage of the model, for there are methods to study the 
determinant of such operators || pj. On the other hand, however, these 
methods can not be directly applyed in the presence of the manifest chiral 
symmetry breaking term in the NJL Lagrangian |l[ [|. We have an even 
more complicated case now, because of isospin breaking. We shall consider 
this problem in the following sections. 
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3 The chiral Ward — Takahashi identities 



One can ask if it is possible to get the symmetry breaking part of the 
bosonized NJL Lagrangian without a direct evaluation of the fermion de- 
terminant in flisp . The answer on this question is positive and hidden in 
Eq. (|l7D . Indeed, all we need to do is just to solve it. Possible solutions 
may be constructed by applying a method which is closely related to the 
one already used in 
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The lack of a direct dependence of chiral trans- 
formations on the space-time coordinates simplifies our task. In particular 
the chiral WT identities can be stated directly by giving the variation of the 
Lagrangian of collective meson fields, £ co ]i, with respect to the infinitesimal 
chiral SU(2) x SU(2) transformations, i.e., 



dL mU — T7 d(T a- 



G 



(19) 



If chiral symmetry is spontaneously broken, the infinitesimal transformations 
are given by Eq.(|IIJ). It implies that the variation of the functional of the 
considered collective meson fields is described by the infinitesimal operator 



5 = 2{a i X l + fcYC, 



(20) 



where the generators Xi and Yi are expressed in the space of meson fields as 
follows 



Xi 



Y 



-ijk 



0-Ku 



a 
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{ai + A,)— + (<7 
dn 



An) A 



& 



7T, 



7T; 



— - 7T — . (21) 



da 



5a i 



The definition of a functional derivative depends on which parameters are 
considered as fixed in the functions, and which are considered as variable. 
Our functions vr a (x) and a a (x) have fixed arguments with respect to global 
chiral transformations, i.e., the variational derivative does not produce a 
delta function of space-time coordinates. In this case the generators satisfy 
the commutation relations 



[Xi,Xj] 
[Xi, Y 3 ] 



^ijkXk , 
^ijkYki 



-ijk 



Xh 



(22) 
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which are the Lie algebra representation of the chiral group. Using the defi- 
nition (PH), we can rewrite Eq . (|19|) as 



-^coll - F i~ -Q e ijk(Tjmk, 

Yi£ C o\\ = Gi = -^(^0^ + rhin ). (23) 

The fact that the WT identities (^) express the variation of a single function, 
£ co jj, gives immediately strong restrictions on the form of this function. In 
addition we have integrability conditions: 

YiGj YjGi ^ijk-^ki 

YiFj — XjGi = eijkGk, (24) 



which tell us if it is possible or not to integrate Eq.(|23|) with the given func- 
tions Fi and Gi. These conditions are fulfilled in our case. It is clear that 
from the WT identities one only obtains general symmetry restrictions on 
the form of the bosonized NJL Lagrangian; in particular, one can use them 
to fix the symmetry breaking part of £ QO \\- Indeed, the first equation in the 
system gives 

where Cs is a chiral symmetric part of the Lagrangian, i.e., XiCs = YiCs = 0, 
and for Q we have 

XiQ = 0, 

YiVL = (rh^i + mi7r + ^-7r a^ . (26) 

Noting that o"o and 7r are cancelled by the action of the operator Xi {X^q = 
XiO-Q = 0), we derive 

G 2GA, 1 ; 



Thus 



£ coll = 2GA^ + ~G~ + U (28) 
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Since we are looking for solutions corresponding to the real vacuum state, the 
linear term ~ cq in the Lagrange density fl28|) must not occur; it would alter 
the Schwinger-Dyson equations. Let us use the freedom inherent in the choice 
of the chiral symmetric piece £$ of £ co n to eliminate the linear scalar field in 
favor of terms quadratic in the meson fields. There is only one chiral invariant 
combination which is bilinear in the meson fields and contains the necessary 
part with a , it is a{a^ + n 2 + 2A cr ). The constant a is then fixed by the 
requirement of cancellation of the linear term in fl28|): a = — rfio(2GAo) _1 . 
As a result, we obtain the solution of chiral WT identities for the bosonized 
NJL Lagrangian in the form: £ co n = Csb + £s, where 

m 3 (n 2 + m (gg + tt 2 ) 
Csb ~ 2GK Z 2GK ' (29) 

We have achieved our stated aim of integrating Eq . (|19|) . Let us turn now 
to the direct calculation of the fermion determinant in Eq.([H|). Contrary to 
the chiral symmetric case, for which the proper-time method yields a correct 
result, we still have to find the way to interpret the fermion determinant for 
the present problem with explicit chiral and isospin symmetry breaking. 

4 Dirac fermion determinant with isospin sym- 
metry breaking 

To evaluate the real part of the fermion determinant we begin with the 
standard proper-time representation (|I|). Although this definition fails when 
chiral symmetry is explicitly broken, it still can be used as a basis for system- 
atic calculations, if one includes the correcting counterterms in the functional 
P. We have also to note that in the present case the standard proper-time 
expansion does not work. There are several reasons for that. First, it is the 
non-commutativity of the quark mass matrix in the heat kernel exponent. 
This property is reflected in resummations inside the standard proper-time 
series. We denote them by T-resummations and show how to organize them. 
Second, this T-resummations do not lead automatically to chiral invariant 
groupings at each order of the P-exponent. One has to search for the invari- 
ants which replace the Seeley-DeWitt coefficients in the asymptotic expan- 
sion of the heat kernel. In practice, this is the most difficult task and implies 
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additional resummations. Third, the gap-equations work against chiral sym- 
metry. One can restore chiral symmetry without afflicting the gap-equations, 
by adding to the Lagrangian a counterterm local in the collective fields. This 
is the general phylosophy we shall adhere to in the remainder of this paper. 
As a result, we obtain the collective Lagrangian, £ co \\, describing the low- 
energy limit of the four quark NJL dynamics, which is in full agreement with 
the WT identities. 

4.1 Proper-time representation for the heat kernel 

The starting point of our analysis in this section is a non-perturbative defi- 
nition of the fermion path integral, or fermion determinant: 

J VqVqexp jz J d A xqDq^ = e iW[n ' a] = det D. (30) 

The functional W[tt, a] = —i log det D is the fermionic effective action and D 
is given by Eq. flT5|) . Of course the fermion determinant defined by Eq. (|30|) is 
only a formal construction. We shall see below that W[7r, er] can be defined 
unambiguously up to a set of local counterterms on the basis of the modified 
Schwinger proper-time representation. We study here only the real part of 
the Dirac fermion determinant. 

Once a definition for det D is chosen satisfying some set of consistency 
requirements, its dependence on the collective meson fields and in particular 
its behavior under chiral transformations can be analyzed. The naive iden- 
tification of the real part of this determinant with the proper-time formula 
Eq.([l|) is problematic for a couple of reasons. First, although the Dirac oper- 
ator D does not include the current quark mass, rh, the transformation law 
of pseudoscalar and scalar fields does. Thus 2 , 

5D E = i[a, D E ) - i{7 5 /5, D E } - 2ry 5 [fn Q f3 + m 3 (p 3 r - 2750^*7*)], (31) 

i.e., the transformation law of the Dirac operator has inhomogeneous terms 
which are proportional to rho and 777,3. As a result we have 

5{D e De) = i[a + 75/?, D E D E ] 

2 For the actual calculations we perform a Wick rotation into Euclidean space-time: 
ixo = £4, ijo = 74, 7 = 7^, {77-, 7s} = — Srs, where r,s = 1,2,3,4. We have chosen all 
7 r -matrices anti-Hermitian: 7^ = — j r . So we get D — > De — i"f r d r — m — a + i757r, where 
75 = 75, d r = (V,<9 4 ). 
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- 4m {/frr + [^e ijk {ai + mJPj + /5 fc vr ] r fe } 

- 4m 3 {75 (a7r 3 - a7rr 3 ) + n (3 3 

+ e 3ik aia k + [P 3 n k + (cr + m )e 3ik ai} r k } . (32) 

It means that 51n|det-D| is not equal to zero and has systematic contri- 
butions proportional to the current quark mass terms m and m 3 . This is 
similar to the case which has already been considered in 0. The explicit 
symmetry breaking terms have to be corrected by the corresponding contri- 
butions from the functional P to lead to a result being in agreement with 
the WT identities. Second, there is a new type of contributions. These ones 
emerge after spontaneous symmetry breakdown and are proportional to the 
isospin symmetry breaking terms ~ (rrid — m u ). Let us discuss these problems 
and their solution in full detail. 

We start by defining the real part of the fermion determinant by the 
formula: 

Re (In det-D) = In | detD| + P. (33) 

The first term here is the Schwinger proper-time regularization for the Dirac 
fermion determinant given by Eq. ([[]). In the case of nonrenormalizable mod- 
els like NJL we have to introduce the cutoff A to render the integrals over 
T convergent. We consider a class of regularization schemes which can be 
incorporated in the expression (P through the kernel p(T,A 2 ). These regu- 
larizations allow to shift in loop momenta. A typical example is the covariant 
Pauli-Villars cutoff JllJ 

p(T, A 2 ) = 1 - (1 + TA 2 )e" TA2 , (34) 

which we use in the present calculations. Our strategy is the following: know- 
ing that the first term in Eq. (^) breakes chiral symmetry in a contradictory 
way, we find the functional P[n, a] by requiring the real part of the fermion 
determinant to transform in accordance with Eq. (|T9|) . The WT identities 
are used for that. The first step in the realization of this program is the 
evaluation of the heat kernel Tr(e~ TR ) with the operator R = m 2 + B, where 

B = —d 2 + ij r (d r a — ij 5 d r n) + a 2 + it 2 + {a, m} + 275(71", a + m). (35) 

We denote by 

m * = Kr + M T3 , K=<±^, M=<^, (36) 
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the square of the mass matrix. The operator R depends on functions of eu- 
clidian coordinates, x r , and derivatives, d r . Therefore, following the abstract 
formalism developed in || , we regard e~ TR as an operator e~ TR acting on a 
fictitious Hilbert space, so that the heat kernel Tr(e~ TR ) reads 

Tr(e" TiJ ) = [ d 4 xti < x\e~ TR \x >, (37) 

where \x > is an eigenvector of a commuting set of Hermitian operators 
x r such that x r \x >= x r \x > and < x\y >= 5(x — y). The Hermitian 
operators, p r = —id r , which are conjugate to x r , obey canonical commutation 
relations: [x r ,j3 s ] = i5 rs , [x r ,x s ] = [p r ,Ps] = 0. Let us take an eigenket \p > 
of p r . Its representative in the Schrodinger representation is < x\p > = 
(27r)~ 2 exp(ix r p r ). Using this plane wave basis one can evaluate the heat 
kernel directly: 

< x\e~ TR \x >= d 4 p < x\e~ TR \p >< p\x > . (38) 



Taking into account the relations: < x'\d r \p >= d' r < x'\p > and < x'\x r \p >= 
x' r < x'\p > one can proceed 

< x\e~ TR \x > = J d 4 p < p\x > e~ TR < x\p > 

P ' e - ipx e- TR e ipx . (39) 



(2^) 4 

We note then that e~ ipx Re tpx — (R — 2ip r d r +p 2 )l, where p r d r = pd, implies 

<x\e- TR \x> = I -^-e- Tp2 e~ T ^- 2ipd h 

(2vr) 4 



" ' I I 'I ■; I i 1 i 

1. (40) 



-e p exp 



— T (k- 2 ^ 



T 



(4tt 2 T) 2 

Finally, combining (|37|) and (^0|), results in the desired relation 

d A pd A x 



Tr(e-™) = / f£^ 2 e-^ +TK h. {exp [-T (Mr 3 + A)} 1} , (41) 
where 

A = B- 2 -^. (42) 

Vt 
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Similarly as in the case of equal quark masses, one has to separate from the 
heat kernal exponent the m 2 piece. However, in the present case m 2 contains 
apart from the isoscalar, K, which we already factorized in Eq.(IT), also the 
non-commuting isovector term proportional to M. This non-commutativity 
is one of the sources of necessary T-resummations in the heat kernel expan- 
sion, which deviate from the standard case and which we consider next. 



4.2 P-exponent and non-commutativity of the mass 
matrix 

We discuss now further the exponent, tr[e~ T ^ Mr3+j4 ^], contained in the heat 
kernel. Here we can not use the standard asymptotic expansion in powers 
of the proper-time T, which would lead finally to the Seeley - DeWitt co- 
efficients. One has first to separate the quark mass part e~™ T3 , otherwise 
it will generate linear terms in o at any power of T. It is not convenient. 
We need a method which leads to the gap-equations already at the first step 
of the asymptotic expansion. This separation of the non-commuting part of 
the squared mass matrix in the heat kernel require the use of the following 
operator identity: 



e 



-T(Mr 3+ A) = e -TMr 3 x p J exp _ J T dsA ( s )\ I ; (43) 



where 

A(s) =e sMT3 Ae- sMT3 . (44) 
The P-exponent above is defined to mean 

P jexp £ dsA(s)^ = 

1 +£("!)" / dsj ds 2 ... ds n A( Sl )A(s 2 )...A(s n ). (45) 
~± Jo Jo Jo 



For our purpose one needs only the first three terms from the right hand side 
of this expression, i.e., 

tr[e- T ( MT3+A )] = tr (e"™^ {1 - TA 
T 2 



+ 4 



A 2 + At 3 At 3 + c{T) (A 2 - At 3 At 3 ) +...}}, (46) 
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where 



c(T) 



2T 2 M 2 V 

1 r TMr 3 
M 2 



2TMt 3 



1 - 2TMr 3 ) 



c?a[M cosh(Ta) — r 3 a sinh(Ta)]. 



(47) 



The details are given in the Appendix. It is clear now that the asymptotic 
heat kernel expansion (|46|) is modified due to systematic resummations gen- 
erated by the terms ~ M. After evaluating the integrals over 4-momentum 
p r one obtains 



Here we took into account that a Lagrangian density is defined up to total 
derivatives. The Y is given by 

Y = ij r (d r a — i^d r 7i) + a 2 + n 2 + {cr, m} + i^^lrr, a + m). (49) 

Finally, the modulus of the fermion determinant may thus be written as 



4.3 Gap-equations against the chiral symmetry 

Several points about Eq.(^) should be clarified right away. The first term 
does not depend on collective fields and has no interest for us. Let us consider 
the second term which is proportional to Y and which we denote as b\. To 
simplify this expression one has to take the integral over T and calculate the 
traces. The last ones include summations over isotopic, color and euclidean 
indices. The integrals over T can be reduced to combinations of some set of 
elementary integrals J n (m 2 ) 




- In | det D E \ = JL | o °° ^ p(T , A 2 )e"™tr ( c -™» {1 — TY 



+ — [Y 2 + Yt 3 Yt 3 + c{T){Y 2 - Yt 3 Yt 3 )\ + 0(Y 3 ) 



T 2 r 



(50) 




n = 0,1,2.... 



(51) 
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These manipulations lead us to the result: 



61 = -^{[M m l) + M m dWl + n 2 a + 2(m a + m 3 a 3 )] 



+ 2[J (m 2 u ) - J (m 2 d )](a a 3 + 7r 7r 3 + m a 3 + m 3 a )} . (52) 



Other terms in the expansion of the integrand in ( |50| ) are ~ Y n , where 
n > 2, and therefore do not contribute to the term linear in a. We consider 
small fluctuations of the system about the asymmetric vacuum state and 
have already shifted the scalar fields correspondingly It means that linear 
terms in the scalar fields should not be present in the Lagrangian. This 
self-consistency requirement can be expressed in terms of the gap-equations: 

m u -m u N C G 2 

-J (mJ, (53) 



m u 2n 2 
m d - m d N C G 



Mm 2 ). (54) 



rrid 2tt z 

Once this has been done, one is confronted with a problem related to the 
fact that the last term in Eq.(p^), linear in the scalar fields, i.e., 



~ ^2l J o( m l) ~ M m d))( m o<73 + m 3 a ) (55) 

vanishes. Let us make clear the essence of the problem. The difference of 
two Jo integrals with different arguments can be written as a sum of two J\ 
multiplied by M and a rest which includes the contributions from J n integrals 
with n > 1 

JoK) - JoK) = -M[JiK) + Mm 2 )] + 0(M 3 ). (56) 

As we soon will see, the first term contributes to the next order of our mod- 
ified proper-time expansion. This fact is very important, because without 
this contribution chiral symmetry would be destroyed at the next order. The 
situation is even more complicated, since the 0(M 3 ) terms will of course 
contribute to all remaining orders of the heat kernel expansion. For the 
same reason relation (|56"D implies also that the term quadratic in the fields 
in Eq . (|52"1) , multiplying the difference of Jo integrals, does not contribute to 
the leading order. Thus in evaluating the fermion determinant on the basis 
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of representation fl5"0|) one has to perform carefully systematic resummations 
in the Schwinger's proper-time expansion. This is an embarrassment but not 
a catastrophy. Indeed, from formula it follows that <5 In | detD E \ ~ rh. 
Therefore the resummations are possible in general, since the total expres- 
sion (0) is already a chiral quasi- invariant, i.e., invariant up to terms ~ rh. 
However, the ansatz introduced by the gap-equations destroys this picture, 
taking away one of the necessary elements. To restore the transformation 
property of the determinant instead of the removed term fl55|) one has to 
add to the functional P a counterterm which is a quadratic polynomial in 
the meson fields, and which has the same transformation property as the re- 
moved expression (|55|). The correction predicted by the WT identities thus 
naturally takes the form 

^[MO - Jo(m 2 d )}-r^[Al(a 2 + 4) + A 2 (<7 2 + vf 2 )]. (57) 
4.4 Extracting quasi-invariants 



Let us turn to Eq.([52]) and consider the first term. This term together with 
the second one in Eq. fll~8| ) gives the part of the meson Lagrangian which is 
responsible for the explicit symmetry breaking. Using the gap-equations one 
obtains 

This result is in obvious contradiction with Eq. (|29|) . But this should not 
surprise us because we already know that the proper-time regularization 
destroys the explicit symmetry breaking pattern of the theory. Fortunately, 
this part of the bosonic Lagrangian can be corrected by the appropriate 
counterterm from the polynomial P which is unambiguously fixed by the 
chiral WT identities. The following circumstances do extremely simplify the 
problem: 

(1) The symmetry breaking part of the NJL Lagrangian must be propor- 
tional to the Jo integrals and not to any other J n . This follows from the form 
of the gap-equations and from the symmetry breaking pattern of the model. 
Thus, all steps in the asymptotic expansion of the heat kernel, excluding the 
first one, must be chiral symmetric. 

(2) It is possible to reduce the problem of finding the explicit form of these 
chiral invariants to the problem of deriving quasi-invariants, i.e., groups of 
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terms which would be chiral invariant if one would add to them the deficient 
terms depending explicitly on the current quark masses. 

The first part of this program requires resummations. The second one 
can be solved by considering corresponding counterterms in the polynomial 
P. 

To see these ideas at work let us restrict ourselves to the easiest case, 
in which the first quasi-invariants will be obtained. We shall carry out the 
necessary arguments in detail only up to the second step in the asymptotic 
expansion; the generalization to the next steps will be more tedious but 
straightforward. Let us consider the third term in Eq. (|50|) , quadratic in Y, 
which we denote by 62- First, one has to perform the T integration. Using 
formula ( pETj ) for c(T), and noting that the part of the integrand ~ sinh(Ta) 
in that expression does not contribute to Eq. ( (50]) , since the trace over isospin 
matrices vanishes, one obtains 



We use now Eq. 

2 



(16vr) : 



rtr 



{[Ji(m 2 u ) + Ji(m; 



Y 2 + Yt.Yt, 



+ 2 



+ M 



JiK) - JiK)]r 3 F 2 
J Q (m 2 d ) - J (m 2 u )] (Y 2 -Yr 3 Yr 3 )}. 



to recast this expression in the form 



where 



+ Q(y 2 -Yt 3 Yt 3 )}, 



(59) 



Y- 



Ji(m 2 u ) 



J^mD] r 3 Y 2 



(60) 



Q 



2M 



K 



Ji(m 2 ) - Ji{m 2 d ) 



M 2 A 4 



K(A 2 + K)- 



+ ... . 



+ 



2MA 4 



(A2 + m2)(A2 + m 2 } 



(61) 



The second term in Eq. (|60D contributes as ~ M and the third one as ~ M 2 . 
It means that these terms contribute to the next and next to the next orders 
of the asymptotic expansion. At the level of the considered approximation 
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one has to take into account only the first term in Eg. fl6"0|) . The trace tr(Y 2 ) 
is equal to 



tr(F 2 



+ 
+ 



AN, 
4 



{d r a a ) 2 + (d r ix a ) 2 + a 2 + 7r 2 + 2(m a Q + m 3 a 3 ) 



n 2 (<Ji + rrii) 2 



(vra + m 3 7r 3 ) 2 + (vrivro + afao + m )) 
m\o 2 + 2m 3 a (7r 7r 3 + a 3 (a + m )) } . (62) 

It is not difficult to recognize a quasi-invariant [a 2 + 7r 2 + 2(m cro + m 3 a"3)] 2 
in the first line. Indeed, the chiral transformations (|13|) leave invariant the 
expression [a 2 + 7r 2 + 2(A o"o + A 3 <r 3 )] 2 which is obtained from the previous 
one by adding the terms proportional to the current quark masses m and 
m 3 . The second and the third lines in Eq.(^) do not have a quasi-invariant 
form. However at this stage there are contributions from b\ (see the term 
~ ((T cr 3 + 7r 7r 3 ) in Eq.Q5^)) and from the polynomial P (see the counterterm 
(|57|)), where in the latter one has to take into account only the parts which 
contribute as [Ji(m 2 ) + Ji(m^)] (see Eq. (|5T)|) . As a result we obtain a quasi- 
invariant structure: 

(tt 2 + a 2 + 2m i a i )(a 2 + n 2 + 2m <y ) - [iri{(Ji + m,) - 7r (cr + m )] 2 . (63) 

We can now adjust this expression to a new one which is invariant under chiral 
transformations by including a corresponding counterterm to the polynomial 
P. The possible counterterm which compensates the non-zero variation of 
( |63|) is easily constructed through the replacement: m a — > A a directly in 
Eq. (|63|) . In agreement with the symmetry requirement we obtain now 

2AiCJi)(cro + 7f 2 + 2A a ) - + \) ~ ttq(^o ' v ' 



(4 + - 2 



A n 



(64) 



instead of (|63|). Finally, if we collect all terms which contribute to the con- 
sidered approximation, we obtain a Lagrange density of the bosonized NJL 
model in the form 

77*0 K + 



coll 



rasK + a 2 ) 
2GA 3 

(d r a a ) 2 + (d r n a ) 2 



+ 



2GA 16tt 
a 2 + tt 2 + 2(A <r + A 3 a 3 



+ 4 



y + - 2 



2A 3 a- 3 )(a- 2 + 7f 2 + 2A a- ) 



(lTi(<Ti + Ai) - 7T ((To + A )) 2 } • 



(65) 
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In closing, we remark that the method developed so far allows for an accu- 
rate derivation of the bosonized NJL Lagrangian which describes the low- 
energy dynamics of Goldstone bosons coupled with scalar fields. This La- 
grangian consistently summarizes the effect of the isospin symmetry break- 
ing, m u ^ m d , and allows for a calculation of its physical consequences. The 
expression fl6"5|) fulfills the chiral symmetry requirements of the fundamental 
NJL quark Lagrangian permitting us to consider it as an effective low-energy 
approximation to the chiral quark dynamics described by the Lagrange den- 
sity (|). 



5 Coupling constants and masses 



For completeness, we would like to discuss here the coupling constants and 
masses of the physical collective modes in the Lagrangian (|65|). We need to 
bring this expression to the standard form, i.e., to diagonalize its bilinear part 
and to introduce the physical states by the corresponding field renormaliza- 
tions. By rewriting the Lagrangian in terms of physical states we complete 
our work in its analytical part. We are not going to give here the numerical 
discussion of the resulting mass formulas for composite mesons or to fix the 
parameters of the model. This numerical part is not relevant for our consid- 
eration here and will be done elsewhere together with the discussion of some 
physical problems. 

Let us consider the part of the Lagrange density fl65|) which is quadratic 
in the fields 



C 



free 



m (cr, 



7T 



2GA, 



2GA r 



{d r a a f + (d r n a ) 



A 



ov°"o 



7T n 



+ A^ 



7T, 



+ 2A A 3 (7ro7r3 + 3aoa3) , (66) 



where 



/ = 



Nc_ 
4tt 2 



Jx{m 2 u ) + Ji(m 



(67) 



There are two kinds of mixing: 00-03 and 7ro-7r 3 . Both of them arise as a result 
of isospin symmetry breaking, and both are diagonalized by the orthogonal 
rotation: 

(<Po\ _ ( cos & sin 6> \ / O \ 
WJ ~ I -sin0 cos 6/J Us/' 



(68) 
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where we have chosen the general notation <fi which means either o or 7r, 
depending on the case under consideration. The rotation angle specified by 
the 0's is equal to 6 S for scalars and to 9 P for pseudoscalars. The quadratic 
form which has to be diagonalized is — (A0q + C0o03 + B(f)f). Then the angle 
of the rotation is fixed by the condition: tan(2#) = C/(B — A). From Eq. fl66|) 
one obtains 

6A A 3 „ 2A A 3 . s 

tan2fl s = — , tan2#„ = — — , 69 

UJ to 



where 



»-*-*+Mt-tY (70) 



For convenience let us rename <pi^ = 0i,2- Having introduced the relevant 
fields, we are now in position to write down the free Lagrangian ( |B"6|) as 

%ee = ^K z 2GA^ " 4 ^ + ^ 

[A s a 2 + B s a 2 + A^l + B p nf) . (71) 
The coefficients with tilde are defined as 

A = Acos 2 + Bsin 2 0-Csm6cos6, (72) 

B = Asin 2 + Bcos 2 6 + Csin6cos6. (73) 

with the following values of coefficients A, B, and C for the scalar and 
pseudoscalar cases respectively (see Eq.(|66|)) 

A p = ^t 3 ~ IA l ^ = ^- /A 3' ^ = 2/A oA 3 . (74) 

The physical states are defined by bringing the kinetic terms in the £f ree to 
the standard form, i.e., 

a a = ga^ h \ n a = g^ h \ Ig 2 = 2. (75) 

Finally these replacements determine the masses to be 



2 

m ^ = ^f:' m l= 2 9 2 Bs, ml o =2g 2 A s . (7(5) 
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Ill 



7Tl,2 



mo9 



ni 



2g 2 B p , 



rn 



""0 



2g 2 A 



(77) 



They are different, for example, from the results obtained in || and R, 12 
They differ also from the ones obtained on the basis of direct calculations 
with Feynman diagrams [13|, [T3| (and references therein). This can be easily 
understood, recalling the general symmetry properties of the model discussed 
in section 2: the chiral transformations of meson fields depend explicitly on 
the current quark masses, see Eq . (|l"3D , but the Dirac operator, ([15]) does 
not. Therefore the Feynman amplitudes calculated directly and solely from 
In det D will never depend on the current quark masses and consequently 
not fulfill the symmetry requirements. To have the correct behavior one 
needs chiral invariant combinations of fields, i.e. the missing current quark 
mass terms have to be included in the correcting polynomial P, via the WT 
identities, as extensively discussed in this paper. 



6 Concluding remarks 

The mesonic degrees of freedom a a , and n a are the elementary collective 
excitations of the non-trivial ground state of the NJL model. The chiral 
dynamics of this collective modes should be consistent with the standard 
picture of a spontaneously broken chiral symmetry, which is naturally linked 
to the explicit symmetry breaking pattern of the underlying four-quark La- 
grangian. In this paper we studied the NJL model with the broken chiral and 
diagonal U(2) flavour symmetry, choosing the quark mass matrix to be in the 
non-degenerate form: m u ^ rhd- We suggested and discussed in considerable 
detail a new method to derive the effective bosonic Lagrangian starting from 
the pure quark action with the chiral four fermion interactions. The collec- 
tive bosonic modes have been introduced through the generating functional. 
It is known that quite often the identification of the path integral over quark 
fields with a determinant is problematic. The case under consideration was 
not an exclusion. The method we developed deals with the accurate defini- 
tion and evaluation of the low energy part of the fermion determinant. We 
studied the real part of the quark effective action and obtained the low en- 
ergy expansion for it. To get the expansion we have chosen the Schwinger's 
proper-time representation of the fermion determinant as a starting point. 
We have shown that in this case one can not use the standard proper-time 
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expansion. The mass matrix of quark fields, which does not commute with 
the mesonic part of the Dirac operator, has first to be factorized from the 
heat kernel. We used the time-ordered exponent to separate the part of the 
heat kernel which depends on the meson fields. In this case the expansion of 
the time-ordered exponent leads directly to the gap-equations. 

On the basis of our recent results |I|, 0] we knew that a proper-time repre- 
sentation failed in the case of manifestly broken chiral symmetry. Therefore 
we used the WT identities to correct the m-dependent part of the fermion 
determinant. In this part of the work we extended the method developed in 
our previous papers to the isospin asymmetric case. However an additional 
new problem arose: we found that the naive expansion of the time-ordered 
exponent has to be modified by additional resummations which touch the 
terms of the effective action which are proportional to the difference of con- 
stituent quark masses — m u . These resummations are necessary to form 
groups of terms which are invariant under the chiral transformations. 

The gap-equation ansatz made the picture even more intricate. We have 
found that because of resummations the contributions of linear terms in 
scalar fields are important for chiral symmetry. As a consequence one has to 
add the correcting counterterm to the effective action to put the ansatz in 
agreement with the WT requirements. 

The method presented in this work is a necessary mathematical tool which 
can be applyed to the more interesting, from the physical point of view, case 
of explicitly broken chiral SU (3) x SU{2>) symmetry. These versions of the 
NJL model have been already studied but without a careful treatment of the 
questions indicated in this paper. 
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Appendix 

Let us obtain here the first four terms in the asymptotic long-wave expansion 
of the operator 

tr [ e - T ( Mr3+A )] . (78) 

This operator is a part of the heat kernel in the expession ([|l]). The result, 
as we already know, has been used in Eq . (f46[) . We begin here by making use 
of the well known formalism of the "time" -ordered exponent, or P exponent: 



tr 



-T(Mt 3 +A) 



tr e- T ^P 



exp 



dsA(s) 



£A, ( 79 ) 



i=0 



where the P-exponent is defined by Eq. fl45|) , and A(s), according to the 
definition (]44|), is equal to 



A(s) 



Ae 



-sMts 



A + T 3 AT 3 + e 2sMT *(A-T 3 Ar 3 



(80) 



One can easily verify that the second term of this expression does not con- 
tribute to the first nontrivial term of the P exponent, for its isotopic trace 
vanishes. Thus the first two terms in Eq.(|46|) are obvious. Let us consider 
now the integral 

T fS 

ds / ds 1 A(s)A(s 1 ) 
o Jo 

ds 



. sM{A 2 + At 3 At 3 + t 3 At 3 A + t 3 A 2 t 3 ) 
47w Jo L 

+ sMe 2sMr3 {A 2 + At 3 At 3 - t 3 At 3 A - t 3 A 2 t 3 ) 
+ 2 sinh(sM)e sMr3 (A 2 - At 3 At 3 ) . 



(81) 



The explicit form for the third term in the expansion of Eq. ( |7D[ ) is given by 



A, 



tr 



2M 



-TMt 3 



ds 



sM(A 2 + At 3 At 3 



+ sinh(sM)e sMT3 (A 2 - Ar 3 Ar 3 )] } 



(82) 



Note that the second term in Eq.(^) does not contribute to A2, by the 
property of the trace, which vanishes. If we integrate over s, we find from 
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( |82"D the third term in the expansion (|46|). The calculation of the coefficients 
Ai becomes a somewhat more and more difficult exercise at every new step 
in (|79"|). To finish this Appendix we give the result of such calculations for 
the coefficient A 3 : 



T 



A 3 = --rtne 



-TMt 3 



^{A,t 3 } 2 + c 2 (T)A{A,t 3 }[t 3j A} 



+ c 3 (T)(A 3 - Ar 3 Ar 3 A)]} , (83) 
where the matricies C2(T) and c 3 (T) may be written as 

c 2 (T) = [l + TMt 3 + (TMt 3 - l)e 2 ™^} , (84) 

c 3 (T) = [e 2 ™^ - 1 - 2TMt 3 - 2T 2 M 2 ] . (85) 
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